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A detailed study of the Inelastic Neutron Scattering (INS) spectra of the high-T c cuprates based 
on the Fermi liquid (FL) picture is given. We focus on the issue of the transformation between 
the commensurate and incommensurate (IC) excitation driven by frequency or temperature. For 
La2-iSr a: Cu04 (LSCO), the condition of small A(0)/vf<i (where a is the lattice constant, and 
henceforth will be set to 1) can simultaneously reproduces the always existing IC peaks in the su- 
perconducting (SC) and normal state, and the always fixed location at temperature or frequency 
change. For Yl^CusOe+i (YBCO), a moderate A(0) /vfcl and proximity of the van Hove singu- 
larity (vHS) at M = (0, 7r) to the Fermi level can reproduce the frequency- and temperature-driven 
shifting IC peaks in the SC state, and the vanishing of the IC peak in the normal state. The 
commensurate peak is found to be more appropriately described as a random phase approximation 
(RPA) effect. We address the conditional peak shifting behavior to a refined consideration on the 
nesting effect which is previously overlook. As a result, both the data on LSCO and the recent data 
on YBCO (on YBa2Cu30e.7 by Arai et al. and YBa2Cua06.85 by Bourges et al.) can be reasonably 
reconciled within a FL picture. We also point out that the one-dimensional-like data by Mook et 
al. on a detwinned and more underdoped sample YBa2Cu306.6 could be due to a gap anisotropy 
effect discussed by Rendell and Carbotte, and we proceed to suggest a way of clarifying it. 
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I. INTRODUCTION 

The aim of this paper is twofold. The first is to give 
a comprehensive study of the essential properties of the 
BCS spin susceptibility, both bare and RPA-corrected. 
The second, basing on the simplest RPA theory we ac- 
count for the present strongly contrasting INS data on 
LSCO § and YBCO §,§. 

Magnetic fluctuation in the high-T c cuprates has been 
long believed to be very intimately related to their su- 
perconductivity mechanism. Therefore various kinds of 
magnetic measurements such as Nuclear Magnetic Res- 
onance (NMR), Nuclear Quadrupole Resonance (NQR), 
and INS have been carried out on the system, and aimed 
for clarifying the connection between the two. Among 
those measurements, INS experiments is the one which 
distinguishes itself from others by its capability to mea- 
sure the fluctuation locally in momentum and frequency 
space and hence making itself an indispensable and im- 
portant tool. Especially recently it provides a great deal 
of new observations, which has yet reached a consensus. 

Essentially there are two dominant features of the ob- 
served INS spectra, the existence of commensurate peak 
at Qaf = (tt, 7r) , and IC peak at Qs = (tt, it ± 5) and 
{it ± S, 7r). The spectra in LSCO j|,fy||] is always IC and 
have a incommensurability independent of lu and T but 
increases with hole-doping. On the other hand, YBCO 



_ and BisSrsCaCuzOg+z (BSCCO) compounds 
: libit both commensurate and IC peaks. Very recently 
a comprehensive INS data on YBCO is obtained j3j§] . In 
the low-temperature SC state, it shows a commensurate 



peak at a particular frequency uj (T — 0). Departing 
from lo (either go below or above uj q ), the peak is split 
up into IC and the incommensurability is continuously in- 
creased. The IC peak is found to appear only above some 
threshold frequency. In the normal state, the excitation 
is always a broad commensurate or weakly IC structure. 
At rising temperature to T Cl the incommensurability of 
the low frequency IC structure at low temperature is con- 
tinuously closed up to commensurate. Such behavior is 
seen in an underdoped S and nearly optimally doped 
YBCO H (while YBCO data on a broad doping range 
are not available). These behaviors are in strong contrast 
to LSCO, and therefore raise an important question as 
whether the INS spectra in YBCO and LSCO are of the 
same origin. 

There are also reports of the observation of the com- 
mensurate resonance in the pseudo-gap phase of under- 
doped YBCO at the same frequency as in the SC phase. 
But that remains controversial ||. At all doping levels, 
its resonance nature is only well-established in the SC 
phase. 

Perhaps the most straightforward and simplest inter- 
pretation of the IC peak is the Fermi surface (FS) nest- 
ing effect (see e.g., Ref. pp|dl^| and references therein). 
Many of the theoretical approaches assume that the sys- 
tem is spatially homogeneous and ultimately end up at a 
FS, could it be a FS of conventional electrons or a FS of 
some exotic fermion such as the spinon. The spin excita- 
tion is constituted by those low energy excitations nearby 
the Fermi surface. Such scenarios are capable of produc- 
ing frequency-shifted IC peaks [IKHI4, and furthermore, 
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the commensurate peak is easily derived either from the 
bare or RPA susceptibility. 

The Stripe ordering picture |l5j is an alternative which 
predicts fixed 1C peaks in spin and charge excitations 
that have incommensurability directly proportional to, 
and depending only on doping. In this picture the system 
is segregated into one-dimensional hole-rich stripes in 
which the holes can freely move, and electron-rich stripes 
in which electrons are ordered antiferromagnetically. The 
electron filling determines the length scale and thus the 
incommensurability. It is believed that the static stripe 
at the 1 / 8-doping [jig , and the 2-to-l-coupled charge and 
spin excitations fll7| are observed. That has boosted a 
lot of works on the "dynamic" stripe which also aims at 
the connection with the superconductivity mechanism. 
A plain Stripe model predicts no commensurate peak. 

Therefore the Stripe picture seems to work well in 
LSCO, while the FL picture works well in YBCO. This is 
intriguing if one adheres to the philosophy that they all 
belong to the cuprate family and therefore should lay in 
a single unified theory. Here we argue that the discrep- 
ancies could be reconciled if the problem is considered 
more carefully. The frequency and temperature driven 
peak shifting present in YBCO but not in LSCO can be 
accounted for by a flatter electronic dispersion near FS 
such that the excitation can be scattered farther from 
it, while electronic dispersion of LSCO is far more steep 
near FS. A refined consideration of the nesting effect, 
which we have termed the dynamic local nesting effect 
is important when the condition A/vp <C 7r is violated. 
In the regime of having shifting IC peak in SC state, the 
IC peak is necessarily vanishes in the normal state. The 
commensurate resonance in the SC state of YBCO and 
BSCCO is addressed as a consequence of the experimen- 
tally observed closeness of the vHS at M to the Fermi 
level ]l8|-p2fl . The effect of the band singularity has been 
shown in many cases to be significant, such as providing 
the vanishing isotope effect, transport anomalies, high 
SC transition temperature, and so on. In short, we have 
argued that the apparently contradictory INS data can 
have a common basis in the FL regime. 

We should also mention the data of Mook et al. 
|H on a detwinned and more underdoped YBCO, 
YBa2Cu306.6- The one-dimensional nature of the data 
where the IC peaks are found only along one of the crys- 
tal axis, is claimed by the authors as a strong evidenc e of 
stripe's existence. Nevertheless we point out in Sec. VII 



that such behaviors can be due to an anisotropy effect in 
the FL picture. The effect is prominent in the frequency 
regime where the data in Ref. was taken. 

This study is also important in the aspect that it pro- 
vides a comprehensive and detailed survey of the major 
consequences from the FL picture, and thus acts as a ref- 
erence to see how much of the experimentally observed 
behaviors are really deviating from it. Comparisons with 
experiments are made where possible. 

Our discussion will be focused on the evolution of the 
peak intensity and location in the parameter space of to 



and T. We organize this paper into sections and a brief 
summary of the findings of each sections or subsections is 
given at their ends where necessary. The first several sec- 
tions are focused on the effect of dispersion vf on the IC 
excitation. Sec. [o| gives the formalism and explains how 
the dispersion at Fermi level vp is modeled. Sec. [II gives 
the bare spectra in SC state, normal state, and at the 
transition. Sec. IV discusses the underlying machinery 
of the formation and behaviors of the IC peaks. Sec. [V] 
studies the effect of inclusion of the Hubbard repulsion 
or AF interaction into the system via RPA. An appropri- 
ate description of LSCO is given. Sec. [v| is devoted to 
the relation between the INS peaks (commensurate and 
incommensurate) and Af-point vHS. An a ppropriate de- 
scription of YBCO is then given. Sec. VII is a discussion 
on the anisotropy effect and th e one -dimensional nature 
in INS spectrum. Finally, Sec. VIII gives a discussion of 
the results and conclusions. 



II. FORMALISM 

What we need in a FL description of the INS spectrum 
is the existence of a well-defined FS and the dispersion 
nearby. All low energy processes are constituted by the 
excitation at the vicinity of the Fermi surface. 

A INS spectrum is proportional to Imx(q, u>) besides 
some Bose-Einstein distribution factor due to the bosonic 
character of the excitation. Here %(q, u>) is the spin sus- 
ceptibility. The BCS bare spin susceptibility is given be- 
low and followed by the RPA-corrected susceptibility. 

The BCS bare spin susceptibility in a one-layer SC 
system (the susceptibility in a normal system could be 
obtained by letting A — ► 0) is 
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where q and ui are respectively the momentum and en- 
ergy transfers. /(-E k ) is the Fermi function and i? k = 
(£ k + A 2 ,) 1 / 2 is the quasiparticle spectrum with £ k and 
A k the band dispersion and SC gap respectively. Hence- 
forth we will write x'o = R- e Xo, Xo = I m Xoi an d likewise 
for x hi Eq. (||). 

We use a two-dimensional tight-binding electronic dis- 
persion 

£ k = — 2i(cos k x + cos k y ) — At' cos k x cos k y — fi, (2) 

where t and t' are the nearest-neighbor (NN) and next- 
nearest-neighbor (NNN) hopping respectively. \i is the 
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chemical potential. In this paper only the d x 2_y2-gap is 
treated (except in Sec. VII) since it is meant to describe 
the high-T c cuprates. It is taken as Ak = A(T)(cosA; x — 
cosk y )/2. The width Y/t in numerical integration is 
taken within 0.002 - 0.008 with slicing 2000 x 2000 or 
4000 x 4000. 

Equation ([j]) describes two kinds of excitation, the 
pair-breaking excitation that excites two quasiparticles 
from the SC condensate and costs energy E^ + -Ek+q, 
and the thermal one-particle excitation that excites a 
quasiparticle from k + q to k which costs energy E^ — 
-Ek+q- The two-particle excitation vanishes in the nor- 
mal state while the one-particle excitation vanishes at 
zero-temperature SC state. 

Here we want to make several important remarks on 
the modeling [by the simple dispersion Eq. (^)] of the 
dispersion gradient at FS (Fermi velocity vf) of a real 
system: 

(i) It is a fact that low energy physics is only rele- 
vant to the part of the dispersion near the Fermi 
level. Therefore though the real dispersion near the 
Fermi level can dramatically deviate from the sim- 
ple tight-binding band (sometimes even shows kink 
in the dispersion p4[), for practical purposes we 
still can model it by choosing an effective scale t 
that gives similar Fermi velocity vf at the Fermi 
level (since vf = \de/dk\Fs oc t). Thus in our con- 
text, A/t is equivalent to A/vf, i.e. a greater A/t 
could mean a smaller vp but does not necessarily 
mean a greater A. Obviously, we do not think "t" 
as the bandwidth of the real dispersion, which is 
not of our concern. 

(ii) On the relative size of A/vp in LSCO and 
YBCO: Though there is a consensus that Aybco ~ 
2Alsco, the uf's are very uncertain |^5|,|4|. We 
therefore take vf as a phenomenological input 
which is to be verified by future ARPES measure- 
ments, and vf is thought as the dominant factor in 
the relative size of A/vf- 

(iii) Note also that since vf (or the effective t) is uncer- 
tain, we will compare all energies (e.g. lo, T, etc) 
only in terms of A. 

The most important correction to the susceptibility 
should be the existence of Coulomb or AF correlations 
between the quasiparticles. Such correlations are be- 
lieved to exist as a residual interaction between the renor- 
malized particles, and they are conventionally treated by 
mean field decoupling as a nontrivial step beyond bare 
theories. Then the susceptibility is written into a RPA 
form as 
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J q is different from U essentially at being positive near 
Qaf and changing its sign to negative near (0,0). The 
RPA-type correction is important because it is perhaps 
the only way to treat correlation effect analytically. 



III. EXCITATION AT BARE LEVEL 
SURVEY OF BASIC PROPERTIES 

We have given in this section a survey of the essential 
properties of the bare excitation spectrum. Their discus- 



sion will be given in Sec. IV 



A. Superconducting State 

The SC state is described by A(T) ^ 0, which results 
in a quasiparticle excitation gap and a SC coherence fac- 
tor. 

Figure [I] shows the frequency-evolution of the spec- 
tra for two typical cases of A/t. We have ignored the 
well-discussed quasielastic node-to-node IC peaks |2(| at 
(tt + 5 0l 7r ± S ) and (n — <5 Q , 7r =b 8 ), and focused on the 
nesting IC peaks at (tt, tt ± S) and (tt ± 8, 7r). In general, 
as lo is increased, the intensity and widths of the peaks 
are increased, and spectrum at Qaf is filled in at some 
frequency near 2A which depends only on band geom- 
etry. When A/t is large, the fill-in is a stepwise jump 
of x"(Qaf,u). This sharp increase of x'o(Qaf,^) has 
an important relation to the resonance at Qaf and will 
be further discussed in Sec. VI. The nesting IC peaks 



are seen to be sharpest at some intermediate frequencies 
before the filling in. 

In contrast to the sharp, well-preserved and not-shifted 
peaks in the A/t <C 1 case [see Fig. |l|(a)], the larger 
A/t calculation [see Fig. 0(b)] shows broader and biased 
peaks that clump into Qaf with increasing frequency. 
Such converging of the IC peak cluster was observed ex- 
perimentally in YBCO and theoretically reproduced by 
us [HH and other groups [ p7|pl]JlO|l . Here we proceed to 



where Vq can be the vertex of Hubbard repulsion U or 
AF interaction 



give its full explanation in Sec. IV. If the saddle vHS 
at k = (0, 7r) is far from the Fermi level, the IC struc- 
ture will be destroyed after the peaks clump in at some 
^/A<2. 

Some stray structures around the IC peaks always exist 
especially in the case of larger A/t and lo/ A. When the 
RPA correction is introduced, they will become relatively 
weaker. 

The main point here is that at lo ~ A, a system with 
moderate A/vp (modeled by a moderate A/t) has its 
nesting peak broad and driven by lo to shift. We will dis- 
cuss in Sec. IV that the shift direction actually depends 
on the orientation of the curvature of the FS near the 
gap nodes. 

B. Normal State 



3 



0.03 



0.02 



0.01 



0.00 



0.09 



0.06 



0.03 



0.00 




(0.5ji,0.5jc) 



(jC,7l) 



(1 .531,71) 



FIG. 1. Frequency evolution of the SC state Xo'(q> w ) with 
different A/t is shown, (a) A/t = 0.03, w/A = 1.0, 1.5, and 
2.0 (from bottom to top); (b) A/t = 0.30, w/A = 0.6,1.0, 
and 1.4 (from bottom to top). In (b), the broad peak shifts 
at changing frequency, and some stray structure away from 
the IC peak is present. Including the RPA-correction [see 
Fig. ^(d)] will make the IC peaks and shifting effect more 
conspicuous. We particularly note that since we mean to 
describe two different systems in (a) and (b), the "A" 
or "t" in (a) and (b) are unrelated quantities. For both 
panels, T = 0, the dispersion is t = 1, t! = —0.25, and 
fi — —0.65. The FS is shown in the inset and the plots are 
along q = (0.5-7T, 0.57r) — (n, n) — (1.57T, n). 

The normal state is described by A(T) = 0, which 
leaves the bandwidth or the Fermi velocity as the only 
energy scale. Eq. (|l|) is then reduced to the Lindhard 
function 



Xo(q,w) = 
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via i?k — ► £k|, £k| = Ck and — £k when £k > and 
£k < respectively. The relation /(— £k) = 1 — /(6t) is 
also used. 

Figure || shows some typical normal state spectra. The 
IC peaks are found to exist only at the regime u>/t <C 1. 
Such destruction of the normal state peaks by frequency 
qualitatively agrees with the experimental fact |28| . In 
that regime, increasing frequency broadens the peak and 
enhances the intensity, but do not shift the peaks. Similar 
to the discussion in previous subsection, u>/t <C 1 in our 
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FIG. 2. Normal state (i.e. A = 0) x"(q, w ) at 
uj/t = 0.04,0.10, and 0.16 (from bottom to top) is shown 
along q = (0.57T, 0.57r) — (n,n) — (1.57T, 7r). Such IC peaks 
exist only in the regime u)/t <C 1. Temperature T = and 
the dispersion is t = 1, t' = —0.25, and [i = —0.65. 

context could simply mean a steep dispersion at FS (i.e. 
large t) and does not necessarily mean a small ui. 

C. Transition between 
Superconducting and Normal States 

In the transition from SC to normal state, the gap 
is continuously diminished and hence the coherence fac- 
tor. The restriction on transition in phase space is also 
relaxed. We will take an empirical relation A(T) = 
A(0)[1 - (T/i;) 4 ] 1 / 2 [see inset in Fig. | at T < T c in 
our discussion. 

If the IC peaks exist in the normal state, at entering 
the SC state [see Fig. |](a)] it w ^ be enhanced when 
w/A(0) > 1 and suppressed if w/A(0) ~ 1, and it will 
also be sharpened for u>/A(0) < 2 since spectral weight 
at Qaf will be gapped out. 

Figure |^(b) shows a remarkable case of transition 
across T c . For a system which was shown previously 
to have frequency-shifted IC peaks in the SC state [see 
Fig. 0(b)], increasing temperature broadens the peaks 
and shifts them towards Qaf- The peaks merge into a 
broad commensurate (or weakly IC structure) structure 
in the normal state. The change is drastic at T < T c 
where A(T) is rapidly changing. 

The essential point here is the following. Given A of a 
system, when vf (model by t) is such that A/vf (mod- 
eled by A/t) is large enough to provide frequency-shifted 
IC peaks, the peaks will also be temperature-shifted and 
vanish in the normal state (since the condition cj/t <C 1 
for the normal state peak to exist is necessarily violated 
at oj ~ A). 



IV. EXCITATION AT BARE LEVEL 
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FIG. 3. Bare spectra Xo(q, u>) with different A(0)/t and 
uj/A[Q) at varying temperature, (a) A(0)/t = 0.03: the 
solid lines show the T — spectra at uj/A(0) = 1.0 and 
2.0 (from bottom to top). The dashed lines has likewise fre- 
quencies but T — T c — 0.25A(0). (b) A(0)/t = 0.30 and 
w/A(0) = 1.0: T/T c = 0,0.80,0.85, and 1.0 (from bottom to 
top). T c = 0.25A(0) and the assumed T-dependence of A(T) 
is shown in the inset. Heating converges the IC peak cluster 
in (b) to a commensurate structure. Note also that neither 
the "A" or "t" in (a) is related to that in (b). The dispersion 
used in both panels is t = 1, t' = —0.25, and /i = —0.65. The 
plots are along q = (0.5-7T, 0.57r) — (n,n) — (1.57r,7r). 

THE DYNAMIC LOCAL NESTING 

The IC peaks are originated from the umklapp and 
inversion symmetry of the dispersion. This section dis- 
cusses the properties of the bare excitation mentioned in 
Sec. [II. Both the uj- and T-driven peak shifting effects 
are essentially due to the dynamic and local nature of the 
nesting effect. 

We start to illustrate the mechanism by a zero tem- 
perature SC system where the bare spectrum is given by 
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The integrand consists of a coherence factor which re- 
flects the non-time-reversal invariance nature of the mag- 
netic measurement, and a S- function that imposes the 
energy conservation rule. 




FIG. 4. Energy contour = uj/2 in the extended BZ 
for u> — 1.0A(0). Some arbitrary local-nesting vectors q/, 
q//, and (\iv are shown. The thick vector between 

above-FS contour piece have better nesting due to the smaller 
curvature. The FS is shown as dotted line, and the dispersion 
is t = 1, t' = —0.25, and /j, = —0.65. The gap here is chosen 
as A/t = 0.50. 



The role of the (^-function should be first discussed 
since it is the most dominant. It limits the contribut- 
ing region of the phase space by the energy conser- 
vation Ek + -Bic+q = uj. At uj ~ 0, the only pos- 
sibility is Ek and E k +q both ~ 0. Therefore q can 
only be vectors connecting the gap nodes, which are 
Qs a = (it + S ,tt ± S ) and (it — S a ,ir ± S a ), with S Q = 



2 sin 



t - i/t 2 - t'fj, I [-2t'\. 



At uj > 0, to satisfy _E k + -E k + q = to, q can be any mo- 
mentum that connects points on the contours _B k = wi 
and -Eic+q = if + ^2 = These are strips of 
area stretched out from the gap node along the FS. For 
most q, it only constitutes a featureless and "incoher- 
ent" background in the q-space. To constitute a salient 
structure out of the background one needs to have extra 
factors that can stack extra weight onto the background. 
A situation that can further provide some weight is when 
u>i = L02 = uj/2 and all the two frequency dependent con- 
tours become £" k = E^+q — uj/2. It spans the whole 
FS when uj — > 2A. Now q can "nest" two contours lo- 
cally by passing through points k Q = (rmr, mt) (where 
m and n are integers) in the BZ. It occurs because the 
w-dependent contour is always locally parallel at k Q + k 
and k G — k due to the dispersion £" k has the k <-> — k 
inversion and umklapp symmetry in the extended BZ. 
This nesting is thus dynamic and local. 

As an example, we consider x"(q, uj) at q G [0 ~ 
27r, ~ 2n] and the integration over k in evaluating x" is 
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FIG. 5. (a) A bare spectrum Xo'CQi^) m the full q-space 
at lu — 1.0A(0). (b) a single branch of "ridge", obtained 
when the integration over k in Eq. (Q) is done only over the 
quadrant (k x ,k y ) G [0 ~ tt, ~ tt]. Four of these 90-degree 
rotated branches superpose to give x"(lj w ) m ( a )- Refer the 
dispersion parameters to Fig. bj(b). 

run over (k x , k y ) G [— tt ~ 71- , —h ~ 7r]. Fig. |] illustrates 
the local-nesting at a frequency lu < 2A. In the inte- 
gration over the first quadrant (k x ,k y ) S [0 ~ tt, ~ tt] , 
those local nesting vectors are shown as q/. Likewise, 
q// for the second, q//j for the third, and q/y for the 
fourth quadrant. 

Figure ||(a) shows the consequence of such local- 
nesting effect. The formation of the structure can be 
illustrated by an illuminating calculation, where the in- 
tegration is done only over the first quadrant (k x ,k y ) G 
[0 ~ tt, ~ %]. Then only one branch of the ridge-like 
structure is obtained [see Fig. ||(b)] and its orientation 
is the same as the contour in quadrant (k x ,k y ) G [tt ~ 
2tt, tt ~ 2tt\. The ends of the ridge are determined by 
Ak = u/2. The locus of the ridge in q-space resembles 
the locus of the contour (which is approximately the FS) 
in k-space but have twice the extension since the nesting 
vectors are twice the vectors describing the locus of the 
FS in k-space [see Fig. |[. 

An IC peak is a result of overlapping of two ridges 
at the symmetrical points Qs = (tt, tt ± S) and (tt ± 
S,tt), where S — 2sin~ 1 [— fi/(2t)], when lu > lus = 



A(T][-fi/(2t)}. For those FS of open topology [see 
Fig. Q, the ridges overlap also at Qs> = (tt + 5' , tt ± 5') 
and (tt — S',tt ± 6') when lu > u>gi = A(T) yj /j,/ 1' , where 
6' = 2 sin yJii/(At'). The above expressions for S ,S, 
and 5' are estimations since the ridges have finite thick- 
ness. 

Frequency wg> = A(T)y/ fi/t' is also the onset fre- 
quency of Xo(Qaf, w), since Qaf also connects equiv- 
alent points of Q«5'/2 (also called the "hot spots"). 

Due to the bending of the ridges [see Fig. ||, one will 
miss the nodal excitation point Q$ o if one scan two nest- 
ing IC peaks at Qa in their nearest distance direction. 
Therefore x" Wlu be seen as always gapped at any q. 
Since this should be more prominent when A/t is small 
(which is likely the case of LSCO), this may explain the 
observation of a momentum independent spin gap by 
Lake et al. on LSCO (§. 

A very peculiar feature which has been overlooked pre- 
viously is the peak shifting behavior, brought by the 
difference of curvature of the contour pieces below and 
above the FS near the gap nodes. This occurs when the 
nodal FS segment has nonzero curvature and the con- 
tour i?k = lo/2 is appreciably opened up from it. In 
the geometry of FS shown in Fig. |^, the above-FS con- 
tour has a smaller curvature and can be nested more 
effectively. That causes the ridges to have biased cross- 
sections, which in turn causes the nesting peaks to have 
their apex biased to some direction [see Fig. |l](b)]. The 
opening up of the contour with lu increases the widths 
of the ridges and peaks, and therefore shifts the apex of 
the broad peaks [see Fig. [|(b)]. The broad and shifting 
behavior of the bare excitation peak are thus necessarily 
coexist. 

Finite temperature in SC states drives the gap and thus 
changes the extensions of ridges both along and away 
from the FS. When T broadens the ridges away from FS, 
it simultaneously shifts the peaks apex like increasing ui. 
Besides, T always causes some smearing in all cases, the 
SC or the normal state. 

When the dispersion at FS (i.e. vf or t) is steep or the 
gap (A) is small, for lu ~ A, the contours are confined on 
the FS because its transverse extension in phase space <5k 
~ A/|de/dk|ps = A/vf <C tt (this is equivalent to the 
condition A/t <C 1). The ridges and IC peaks formed 
are slim and no shift of peak is possible. For the peak 
to shift, a finite curvature of the FS is also necessary to 
provide the difference of curvature between the above- 
and below-FS contours. 

The coherence factor has values between and 2, and 
that can provide modification of the intensity as pointed 
out in the literature. The condition £k — > at the 
FS tends to maximize the spin coherence factor to 2 at 
AkAk+q < 0, and minimize it to at AkAk+ q > 0. As 
our nesting vectors has AkAk+ q > 0, the ends of ridges 
are smoothly suppressed by the cutoff Ak = lo/2. In 
the case of larger A/t, the condition £k — * is not well 
satisfied at finite to/ A since the contour is opened up 
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appreciably from the FS. This weakens the suppression 
effect before cutoff, and the ridges will overlap at more 
sizable values and gives a higher contrast of the nest- 
ing IC peaks. It also lowers the threshold cj/A for the 
emergence of the peaks. 

The normal state is the lift of the energy restriction 
by the gap, and the disappearance of coherence factor. 
The ridges now always span the whole q-space and the 
nesting peaks always exist in principle. But due to re- 
placing a stronger energy constraint (due to two-particle 
excitation) by a looser constraint (due to one-particle ex- 
citation), the transition now floods over regions far away 
from the FS and makes the widths of ridges dispersed 
and the nesting peaks less distinguishable. Especially 
now the transition near Qaf is filled high for hole-like 
FS. 

Albeit the underlying cause of formation of the ridges 
and nesting peaks (umklapp and k-inversion symmetry) 
is irrelevant to the existence or not of the SC gap, there 
are still some upper limits on uj/A or uj/t for the IC 
peaks to exist. Since increasing tu/t broadens the ridges, 
the ridge and peak structure will be always destroyed at 
very high toft. The same is true for very high T/t where 
the excitation is deconfined from the FS. It is easier for 
the peak to survive in the SC state because the near Qaf 
transition is kept away at u/A < 2 [see Fig. 0(b)] . 

To summarize, at low uj/t and T/t, the nesting IC 
peaks should always exist in bands with umklapp and 
inversion (k «-> — k) symmetry. In the SC state, it can 
be driven by w or T to shift when A/t is appreciable and 
the nodal FS is curved. 



V. RPA-CORRECTED EXCITATION 

In this section we will discuss the effect of RPA cor- 
rection on the frequency-driven shifting effect of the IC 
peaks. We start with a brief discussion of the essential 
features of the RPA spectrum. 

The RPA spectrum is given by 

r, ( \ = Xo(q^) , 7 s 

X ^ ' [1 - ^(q, u)Y + [^(q, u)} 2 ■ 1 ' 

It can be viewed as a scaling of x'o by the Stoner factor 
[(1 — Vqx'of + (V q x") 2 ]i which is greater or equal to zero. 

Figure @ shows that Stoner factor s = [(1 — Vx' ) 2 + 
(Vx") 2 ] versus Vx'o f° r various Xo/x'o- ^ summarizes all 
features of the RPA correction in a lucid way. For in- 
stance, at the regime of small x'Ux'oi RPA spectra with 
V < 1 /x'o is just a scaling (i.e. s ^ 0, could be magni- 
fication or suppression) of the bare spectra; while RPA 
spectra with V ~ 1/Xo are resonances (i.e. s ~ 0) with 
damping x'o- For large Xo/x'oi s i s greater than 1 in most 
regions, which means that for most choices of V the RPA 
correction has only suppression effect. This can happen 
in systems with large x'o due to band singularity. 




FIG. 6. The Stoner factor [(1 - Vx'o) 2 + (Vx'o) 2 } is plotted 
against Vx'o for x'o /x'o = 0, 0.5, 1.0, 2.0, and 10 (from bottom 
to top). 

In reality, the RPA resonance could be unlikely to 
occur in systems with well-behaved properties near the 
Fermi level by the following reason. For any stable sys- 
tem, Vq should be sufficiently away from the instability 
boundary Vq — l/Xo(q, w = 0) [see Fig. [j] for an exam- 
ple]. When there is no singular sources such as the van 
Hove singularity near the Fermi level, x"(q, to) changes 
smoothly with u> and x' Q (q, cj) is inert to change of u> 
(due to the Kramers-Kronig relation). Theoretically a 
fine tuning of the interaction strength is then required 
to get a resonance. In such systems, absence of instabil- 
ity at zero to is likely to imply absence of resonance at 
nonzero u). 

Figure || illustrates the effect of RPA correction on the 
peak shifting behavior. We have considered antiferro- 
magnetic interactions J q which do not cause resonance 
(weak coupling). In the case of moderate A/t, sharp lo- 
cal minima of 1/Xo [Fig- !( c )] appear at the sharp edges 
of the bare IC peaks due to the Kramers-Kronig relation. 
This is nontrivial as it results in more spectacular shifting 
of the IC peaks [compare Fig. ||(d) and Fig. 0(b)]. On 
the other hand, in the case of small A/vp, the fixed in- 
commensuration behavior is maintained [Fig. ^|(b)] due 
to the "quiescent" l/x' [Fig- ||( a )]- We conclude that 
the weak coupling RPA correction do not alter the orig- 
inal behaviors of the IC peaks, and the LSCO data can 
be described by small A/vp and weak coupling RPA cor- 
rection. 



_VI. EXCITATION SPECTRUM AND 
A/-POINT VAN HOVE SINGULARITY 

This is a RPA spin excitation theory for the recently 
reported INS data on YBCO. For simplicity, we have 
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FIG. 7. A typical SC state -l/x' (q,w = 0) in 
q-space. The no-instability bound of the interaction 
strength is < l/Xo(q,w = 0). For V q = J q , 

1^1 < l/Xo(q ~ Qaf,u = 0) ~ 2.2t in this case. The plot has 
T = 0, A/t = 0.03, and band dispersion t = 1, t' = -0.25, 
and = —0.65. 

mimic the realistic extended-saddle- vHS at the proximity 
of Fermi level |ll|^2| by simple saddle-vHS just beneath 
Fermi level. It is believed that in RPA theories they will 
be qualitatively the same. Nevertheless its extended na- 
ture is crucial for bare theory to have quantitative agree- 
ment with experiments p7| ]. 

A number of comprehensive theories including both 
commensurate and IC excitation [ p^pl]j27|Jl2] -p^| were 
explored since the discovery of the behavior in YBCO. It 
is important to note that all of them arc Fermi- liquid-like 
theories, and most are RPA theories Jl^,|ll|,|l3|,y) . The 
possibility of a bare theory was investigated by Abrikosov 
p7| , p9t and us [Q. The commensurate peaks in bare 
theories are generally weaker, broader, and occurs at 
uj > 2A. That is in contrast to the much stronger 
and sharper one in \" that occurs at uo < 2A. How- 
ever, recently it has been shown that incorporating the 
extended-vHS |27j into bare theory can also give peak 
intensity and width in reasonable agreement with exper- 
iments. Attempting to rely on the criterion uj being less 
or greater than 2A for choosing a theory also encounters 
the subtleties in the determination of the SC gap. Es- 
pecially in the underdoped samples where there are still 
controversies on whether the SC gap and the pseudogap 
are the same gap |30|-|32|| . Therefore either lo > 2 A (see 
discussion in Ref. |29|) or uj < 2A (e.g. in the under- 
doped compounds) can be occasionally favored. As a re- 
sult we believe that at present there is no clear evidence 
here to favor either scenarios. Recently the IC peaks 
have received more attention because of their frequency- 
shifted behavior. Most of the present scenarios for it are 
based on the nesting effect which depends mainly on the 
nodal FS, except one recently proposed by Abrikosov |27]] 
depends on the band mass at the extended-vHS region. 
The discussion of the temperature-shifted IC peaks has 
been missing from the literature. 
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FIG. 8. RPA effect (with no resonance) on two systems 
with small and moderate A/t. All graphs except (a) are plot- 
ted along q = (0.57T, 0.57r) — (7r,7r) — (1.57r,7r). For the case 
of A/t — 0.03 and \ J\/t — 1.0: (a) l/Xo(ci, w ) is shown [along 
q = (0.757r,0.757r) - (7r,7r) - (1.257T, tt)] at uj/A = (dashed 
line), 1.0, 2.0 and 1.5 (solid lines at Qaf from top to bottom); 
Correspondingly, (b) x' (c±, w) is shown at uj/A — 1.0, 1.5 and 
2.0 (from bottom to top). For the case of A/t = 0.30 and 
\J\/t = 1.0: (c) l/x' (q,w) is shown at uj / A = (dashed 
line), 0.60, 1.0 and 1.4 (solid lines from top to bottom); Cor- 
respondingly, (d) x"(q, w) is shown at uj/A = 0.6, 1.0 and 1.4 
(from bottom to top). Note the sharp minima of 1/Xo in ( c ) 
and its discussion in the text. Temperature T = and the 
corresponding bare spectra were shown previously in Fig. [I]. 
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With a sufficiently strong correction strength, it is al- 
ways possible to have a RPA resonance at Qaf- The rea- 
son is since any susceptibility function or interaction ver- 
tex on the lattice is symmetric at Qaf + <5q — > Qaf — 
they always have an extremum at Qaf [e.g. see Fig. |?]] 
and thus one can always have 1 — Vx' — locally at 
Qaf- 

The dispersion taken here to mimic the YBCO dis- 
persion has a hole-like FS and a simple saddle-vHS at 
M = (0, 7r) slightly beneath the Fermi level (£ w= jj/A = 
—0.5). The interaction is chosen as an antiferromagnetic 
J q with \ J\/t = 1.6. At T = 0, a commensurate reso- 
nance occurs at lu/ A = 1.5. The momentum dependence 
of the RPA spectrum and 1/x'o are shown in Fig. [|, at 
different frequencies and zero temperature. The onset of 
Xo(QaFj w ) an d the minimum of 1/x' (QaFi w) occurs 
at lu/A = 1.77. 

A new feature appears at u> > lu d (T = 0) due to the 
shallow saddle- vHS. The discommensuration is seen to 
increase at increasing frequency [see Fig. ^|(d)]. It is also 
due to a dynamic local nesting effect similar to that at 
lu — > lu from below. Owing to the shallow flat band 
at M, the nesting contours = lu/2 at high frequen- 
cies necessarily have smaller curvature at those segments 
below FS, and can be more effectively nested locally [see 
Fig. [To) . The peak profile here is weak because the below- 
and above-FS contours are too far apart to coherently 
form a ridge. The RPA correction has an effect similar 
to that at lu < cu (T = 0). 

Replacing the saddle-vHS by an extended-saddle-vHS 
do not affect the IC excitation at lu < lu q since only the 
nodal FS is of concern. Neither the diverging discom- 
mensuration at lu > lu will be affected since the switch- 
over of the smaller and larger curvature energy contours 
should also occur. 

The large change in 1/x'o at a;/ A = ~ 2 is a re- 
sult of the large jump at the onset of x"{Qaf (at 
lus'/A(T) = \] \ijt' ~ 1.8). The onset is sharp and large 
due to the enhancement by the spin coherence factor (at 
AkAi c +Q 4F < 0), and the picking up of states around the 
shallow vHS at M (at a depth of £, k= jj — 4t' — /i) by the 
S- function [see Eq. (g)] . This widens the window and low- 
ers the threshold for the interaction strength giving res- 
onance. This may partially account for the existence of 
the commensurate peak in YBCO and BSCCO, but not 
in LSCO. The RPA effect at Qaf is effectively turned off 
when lu is increased beyond the maximum of x" (Qaf, w), 
since 1/x' {Qaf,u) recedes rapidly from zero then [see 
Fig. ^(c)]. Since the onset frequency scales with A(T), 
the RPA commensurate resonance is necessarily softened 
by temperature. But the softening could be very little. 
The reason is that as A(T) gets smaller, divergence in 
Xo(QaF)W) gets smaller and will fail to fulfill the res- 
onance condition x' (Qaf,^) = 1/V(Qaf) for a given 
V(Qaf)- The peak could cease to exist before its soft- 
ening is really detected. The case for the softening of the 
commensurate peak in a bare theory is similar [p.2[ . 




(0.7571,0.7571) (7C.7C) (1.2531,31) 




FIG. 9. Frequency-driven peak shifting at different fre- 
quency regimes : (a) l/% (q, a>) at ui/A — 1.0,1.2, and 
1.5 [solid lines from top to bottom]; (b) x"(q,lv) at 
lu/ A — 1.0, 1.2, and 1.5 [at decreasing incommensurability 
along (7r,7r) - (1.257r,7r)]. (c) l/x' (q,w) at lu/A = 1.5,2.3, 
and 2.5 [solid lines at increasing "incommensurability" 
of the minima along (n,iv) — (1.5tt, 7r)]; (d) x'io.^) a t 
lu/A — 1.5,2.3, and 2.5 [at increasing incommensurability 
along (n, n) — (1.57T, 7r)]. The dashed lines from top to bottom 
in (a) and (c) show 1/x'o a t lu/A — and 1.77 respectively, x" 
in (b) and (d) are calculated with \ J\/t = 1.6, and the com- 
mensurate resonance occurs at lu = lu (T = 0) = 1.5A(0). We 
have truncated the strong commensurate peak for obvious rea- 
son. T = 0, the quasiparticle dispersion is t = 1, t' = —0.20, 
H = -0.65, and A = 0.30 (gives £ k=I7 /A = -0.5). The inset 
shows the FS. The broadening taken for all the calculations 
in this section is F/t — 0.008. 
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FIG. 10. Energy con- 

tour E k = uj/2 for uj = 2.5A(0) > uj (T = 0). The thick 
local-nesting vector between below-FS contour has more pro- 
nounced nesting effect (compared to the thin vector) due to 
the smaller curvature. We remind that the dispersion here has 
a saddle-vHS near to Fermi level (refer the exact parameters 
to Fig. ^|). The FS is indicated as dotted line. 

Figure |ll| plots the temperature-evolution of a low fre- 
quency spectrum in momentum space. The most no- 
table should be the softened commensurate resonance at 
uj < u> (T = 0). When the system is cooled down from 
above T Cl the broad and weak commensurate peak in the 
normal state grows into a sharp and strong peak at some 
temperature just below T c , then it is split up into IC 
peaks and recede from Qafi to a fixed location at low 
temperature. The weak IC structure at u> > uj (T = 0) 
in SC state is found to be simply smeared off to a broad 
commensurate at warming up to the normal state p3| . 

Note that the observed relative excitation intensity in 
the course of the evolution [see Fig. 4 in Ref. g]] can only 
be obtained correctly in x" (Qj <*0 via an appreciable in- 
teraction strength. Though the "shoot and split" behav- 
ior of the uj < u> (T = 0) spectrum at cooling down is 
already contained in a bare theory invoking vHS effect 
for the commensurate peak |T3], the normal state com- 
mensurate structure always incorrectly has an intensity 
higher than the low temperature IC structure in x'o (l> ^-0 
[see, for example, a similar case in Fig. |](b)]. 

A weak and broad commensurate structure is always 
seen in the normal state at all frequencies [see Fig. |lj] 
as what is observed in experiments j^J^] . It is due to the 
enhancement by the saddle- vHS and RPA correction. For 
the model system considered here, its intensity reaches 
a maximum at some w/A(0) ~ 0.4 as a result of the 
competition between x'o an( i RPA enhancement. 

In summary, the diverging incommensurability at u> > 
lj in the SC state is accounted by the approaching of 




(0.757t,0. 75tc) 



(1.2571,71) 



FIG. 11. Temperature-driven peak shifting at low fre- 
quency regime : At u) = 1.0A(0) < uj {T = 0), (a) l/x' (q,w) 
is shown at T/T c = 0, 0.8, 0.88 (solid lines from top to 
bottom), and 1.0 (dashed line); (b) x"( t lj a; ) i s shown at 
T/T c = 0,0.8,0.88 (solid lines at decreasing incommensura- 
bility), and 1.0 (dashed line). The weak commensurate struc- 
ture in the normal state is seen to grow into a strong peak at 
T = 0.88T C < T c , then splits into IC peaks which recede from 
Qaf at T — > 0. The strong commensurate peak is truncated 
for obvious reason. T c = 0.25A(0), and the T-dependence of 
A(T) is assumed as the same as in Fig. gb). Other parame- 
ters are referred to Fig. ^. 

M -point vHS to the Fermi level; the normal state broad 
commensurate structure at any frequency also finds a 
natural explanation in it. The recent observation on the 
relative excitation intensity in the normal and SC state 
prefers a RPA theory in the resonance regime more than 
the bare theories, though the gross features of the bare 
and RPA spectra are the similar. Including an under- 
doped sample (by Arai et al. 0) and a nearly optimally 
doped sample (by Bourges et al. [pi), the converging and 
diverging shifting of the IC peaks has a transparent ex- 
planation in the FL picture. 



VII. ANISOTROPY EFFECT 

Anisotropy effect in the high-T c cuprates was discussed 
before by Rendell and Carbotte J34[. The effect is be- 
lieved to be important in systems such as the orthorhom- 
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FIG. 12. The broad commensurate peak of normal state 
x"(q,w) at w/A(0) = 0.2,0.4 (dashed lines from bottom 
to top), 0.8,1.4, and 2.0 (solid lines from top to bottom). 
T = T c — 0.25A(0). Refer those not stated parameters to 
Fig- § 

bic crystals or the chained- YBCO. In those systems, 
anisotropy may either exist in the band dispersion or en- 
ter into the SC order parameter, or exist in both. It was 
pointed out in Ref. j34| that if only one of those two kinds 
of anisotropy exist in the system, it can be distinguished 
by the criterion that the gap anisotropy exist only in the 
SC state while the dispersion anisotropy persists into the 
normal state. Some important effects of it which can 
be seen in INS spectrum are, if there is an admixture 
of s-component in the d-wave gap then there should be 
a rendering of the four-fold symmetry of the quasielas- 
tic node-to-node IC peaks to two-fold, and regardless of 
the anisotropy source a difference of the nesting IC peak 
intensities on different crystal axis should be observed. 

We will discuss here a case of anisotropy in the gap, 
the existence of an s-component and the SC order pa- 
rameter has the form Ak = A(T) [ask + (1 — a)dk] , where 
Sk = 1 and dk = (cos k x — cosk y )/2. The existence of an 
s-component is suggested by the tunneling experiments 
p5| [37| and its percentage is at a ~ 0.2. It was discussed 
in Ref. |Q at low frequencies, but in order to clearly dis- 
tinguish the effect we think that it is necessary to study 
its ^-dependence. 

Figure |l3| shows the contrast of the nesting IC peaks 
along different crystal axis at different to, on a system 
with reasonable size of s-component in the gap. The 
difference in peak intensities is seen to be most promi- 
nent at low to, and rapidly decreased at increasing fre- 
quency. Since the data by Mook el al. (on a detwinned 
orthorhombic YBa2Cus06.6) |23| is at a somewhat low 
frequency 24 meV ~ A(0), we believe data on a wider to 
range is desirable in order to make definite statement on 
the origin of the one-dimensional nature. 

The difference in peak intensities can be easily ex- 
plained by the picture given in Sec. IV. The twofold sym- 
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13. Anisotropic 
to) (a) in the full q-space at w/A 



FIG 
X"(q 

q = (-7T, 0.571") — (-7T, 7r) 

(from bottom to top) 



SC state 
1.0, and (b) along 
- (1.5tt, it) at lu/A = 0.6,0.9, and 1.2 
The difference in height of the peaks 



on different axis is more spectacular at low frequencies. The 
SC order parameter is taken as Ak = A[0.2sk +0.8dk], where 
A = 0.10, s k = 1, and <f k = (cosk x — cosfe H )/2. T = and 
the dispersion is t = 1, t' = —0.25, and /i = —0.65. 

metric spectrum has the node-to-node excitation closer 
to one of the two crystal axis, and excitation away from 
the nodal excitation is suppressed by the coherence fac- 
tor, therefore the peaks on one of the axis has lower in- 
tensity. This can only dominate at frequencies as low 
as ~ A and excitation along the FS has not been fully 
opened up. 

We have presented here an anisotropy effect suggested 
by the tunneling experiment and it can introduce one- 
dimensional feature into the INS spectrum (other possi- 
bility is also suggested in Ref. ||]). The frequency de- 
pendence of the contrast of IC peak intensity should act 
as a criterion to compare with the Stripe interpretation 
which has no frequency dependence. Furthermore such 
gap anisotropy effect should vanish in the normal state. 
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VIII. DISCUSSION AND CONCLUSIONS 

We have presented a thorough study on the basic prop- 
erties of the spin susceptibility and have summarized our 
findings and comparisons with experiments at the end of 
each section. 

Some possible departure of our conclusions from real- 
ity should be mentioned. The presence of impurity in 
real systems destroys the exact translational symmetry 
and mixes states of different momentum, which are also 
the energy eigenstates. Such process necessarily relaxes 
the energy conservation selection of the transition region 
in phase space. The effect should be significant in low 
energy processes in the SC state where there is inhomo- 
geneous gapping out of phase space by the d-wave gap. 
For example, it may pull down the lower threshold for the 
existence of the nesing IC peaks |Q . We have also taken 
a simple d-wave in the discussion which is just meant to 
describe a symmetry with four nodes and change of sign. 
The actual k-dependence of the gap should exhibit some 
deviation from pure d-wave since other factors such as 
interlayer tunneling []39f or correlation effects etc. may in- 
fluence. The importa nt ca se of possessing a s-component 
was discussed in Sec. VII. The realistic dispersion at the 
FS also shows deviation from our simple tight-binding 
band, the most well known of all should be the presence 
of extended- vHS. Therefore all numerals we give are qual- 
itative. 

To a qualitative level, the "inconsistency" between INS 
spectra of LSCO and YBCO as mentioned in Sec. | can 
find a consistent picture in the FL interpretation. They 
could be systems in different regimes of the ratio A/vf. 
The observations on LSCO [as discussed in Sec. |l| and 
Sec. [vj: the wide existing IC excitation with incommen- 
surability depends only on doping, can be ascribed to a 
small A/vp of the system. The observations on YBCO 
[as discussed in Sec. [II and Sec. VI \: including the u>- 
and T-driven shifted IC peaks in SC state, and normal 
state broad commensurate peak, can be all ascribed to 
an appreciable A/vp of the system. The diverging in- 
commensurability at to > lo in SC state are ascribed to 
the approaching of the M-point vHS to the Fermi level. 
Though bare and RPA theories could be qualitatively 
similar, RPA theories can describe the temperature evo- 
lution of the excitation intensity more adequately. A pe- 
culiar but straightforward conclusion for the commensu- 
rate resonance is it can occur at u> < u) (T — 0) when 
T < T c , either in the bare U) or RPA theories. This nat- 
urally explains the observed enhancement of scattering 
intensity at Qaf in that regime [see Fig. 4(c) in Ref. M]. 

We would like to point out that while the observed 
commensurate peak in the SC state could easily find an 
interpretation as a resonance, as it occurs at a small fre- 
quency window; the lower intensity broad commensurate 
peak in the normal state which is observed at any fre- 
quency is hard to be interpreted in terms of a resonance. 
In our scenario this is a natural consequence of the over- 



lap of broad IC peaks, with further enhancement by the 
M -point vHS and RPA-correction. 

On the ratio A/vp of LSCO and YBCO, there should 
be a consensus that the SC gap of LSCO (~10 meV) is 
a fraction of that of YBCO (~30 meV). Due to the low 
resolution in the determination of electronic dispersion 
(the resolution in ARPES measurements is about 10-40 
meV), presently vf in most material could not be pre- 
cisely determined po| , p5|j24| and neither the ratio A/vf- 
Therefore in this paper we have meant only a first-step 
qualitative study, and did not predict any precise figure 
for vp of different systems. But it is still important to 
note that in our scenario, vf of YBCO should be at least 
several times smaller than that of LSCO . 

In the literature, there are also other approaches 
0,[ll]|27]|l| to the issue. As we do, Ref. @ also simul- 
taneously accounts for both the LSCO and YBCO data, 
but does not elaborate on how the difference in the be- 
havior of the IC peak arises. Furthermore, it is predicted 
that a commensurate peak should also exist in LSCO. 
Ref. Q addresses only the YBCO data by using a re- 
fined tight-binding band. It is also not shown explicitly 
that how the IC peaks are shifted. Moreover, the relation 



to the LSCO data is not mentioned. Ref. 27 obtains the 
commensurate and shifting IC peaks from the quasiparti- 
cle masses at the extended-saddle-vHS region, while the 
data of LSCO is also not mentioned. In our present work, 
we assume that the IC peaks in LSCO and YBCO have 
the same origin, and have related the different behaviors 
of them to their individual Fermi velocities vp- An ex- 
plicit explanation of the shifting behavior is also given. 
The approach in Ref. [[i"lf is very different from ours as 
it interprets the shifting of the IC peak as the downward 
dispersion of a spin wave mode. Our approach can also 
easily address the shifting of IC peak caused by temper- 
ature change, which is so far only discussed by us. 

There are some bewildering stems from our work. In 
the underdoped cuprates, anomalies such as the pseudo- 
gap phenomenon or the destruction of FS [|42"f are seen 
near the maximum gap region and it is likely that the FL 
picture is to breakdown at that region. But now the FL 
behavior near the maximum gap region is an important 
ingredient in our scenario on the YBCO, and our cover- 
age includes both the underdoped and optimally doped 
YBCO. It seems that to some extent the FL behavior is 
still retained in the response to INS measurements. 

In conclusion, albeit the FL picture on the cuprates is 
recently quested by many experiments, most behaviors 
of the recent INS data are still well described in terms of 
the traditional Fermi liquid picture. 



ACKNOWLEDGMENTS 

KKV acknowledges the scholarship from NSC of Tai- 
wan under grant No.89-2112-M-003-009. We thank P. 
Bourges et al. for sending us their results prior to publi- 



12 



cation. We also want to thank A. Abrikosov for sending 
us his preprint as well as useful correspondences. 



K. Yamada, C.H. Lee, K. Kurahashi, J. Wada, S. Waki- 
moto, S. Ueki, H. Kimura, Y. Endoh, S. Hosoya, G. Shi- 
rane, R.J. Birgeneau, M. Greven, M.A. Kastner, and Y.J. 
Kim, Phys. Rev. B 57, 6165 (1998). 
M. Arai, T. Nishijima, Y. Endoh, T.Egami, S. Tajima, 
K. Tomimoto, Y. Shiohara, M. Takahashi, A. Garrett, 
and S.M. Bennington, Phys. Rev. Lett. 83, 608 (1999). 
P. Bourges, Y. Sidis, H.F. Fong, L.P. Regnault, J. Bossy, 

A. Ivanov, and B. Keimer, Science 288, 1234 (2000). 

B. Lake, G.Aeppli, T.E. Mason, A. Schroeder, D.F. Mc- 
Morrow, K. Lefmann, M. Isshiki, M. Nohara, H. Takagi, 
and S.M. Hayden, Nature 400, 43 (1999). 

J.M. Tranquada, J.D. Axe, N. Ichikawa, A.R. Mooden- 
baugh, Y. Nakamura, and S. Uchida, Phys. Rev. Lett. 
78, 338 (1997). 

P. Dai, H.A. Mook, and F. Dogan, Phys. Rev. Lett. 80, 
1738 (1998). 

H.F. Fong, P. Bourges, Y. Sidis, L.P. Regnault, A. 
Ivanov, G.D. Gu, N. Koshizuka, and B. Keimer, Nature 
398, 588 (1999). 



H.A. Mook, F. Dogan, and B.C. Chakoumakos, cond 



mat/9811100 



P. Bourges, B . Keimer, L.P. Regnault, and Y. Sidis, cond- 



mat/0006085 



Y.-J. Kao, Q. Si, and K. Levin, Phys. Rev. B. 61, R11898 
(2000). 

M. Norman, Phys. Rev. B 61, 14751 (2000). 



K.-K. Voo and W.C. Wu, |cond-mat/9911321 



F. Onufrieva and P. Pfeuty, :ond-mat/9903097 



J. Brinckmann and P. Lee, Phys. Rev. Lett. 82, 2915 
(1999). 

V. Emery and S. Kivelson, Physica C 235-240, 189 
(1994). 

J.M. Tranquada, B.J. Sternlieb, J.D. Axe, Y. Nakamura, 

S. Uchida, Nature 375, 561 (1995). 

H. Mook and F. Dogan, Nature 401, 145 (1999). 

J. Ma, C. Quitmann, R.J. Kelly, P. Almeras, H. Berger, 

G. Margaritondo, and M. Onellion, Phys. Rev. B 51, 
3832 (1995). 

K. Gofron, J.C. Campuzano, A. A. Abrikosov, M. Lin- 
droos, A. Bansil, H. Ding, D. Koelling, and B. Dabrowski, 
Phys. Rev. Lett. 73, 3302 (1994). 

D.M. King, Z.-X. Shen, D.S. Dessau, D.S. Marshall, C.H. 
Park, W.E. Spicer, J.L. Peng, ZY. Li, and R.L. Greene, 
Phys. Rev. Lett. 73, 3298 (1994). 

D.S. Dessau, Z.X. Shen, D.M. King, D.S. Marshall, L.W. 
Lombardo, P.H. Dickinson, A.G. Loeser, J. DiCarlo, C- 

H. Park, A. Kapitulnik, and W.E. Spicer, Phys. Rev. 
Lett. 71, 2781 (1993). 

K. Gofron, J.C. Campuzano, H. Ding, C. Gu, R. Liu, B. 
Dabrowski, B.W. Veal, W. Cramer, and G. Jennings, J. 
Phys. Chem. Solids 54, 1193 (1993). 



[23] 
[24] 

[25] 



[26] 
[27] 
[28] 

[29] 
[30] 

[31] 

[32] 



[33] 
[34] 
[35] 



[36] 

[37] 
[38] 
[39] 
[40] 



[41] 
p] 



H.A. Mook, P. Dai, F. Dogan, and R.D. Hunt, Nature 
404, 729 (2000). 

P.V. Bogdanov, A. Lanzara, S.A. Kellar, X.J. Zhou, E.D. 
Lu, W.J. Zheng, G. Gu, J. -I. Shimoyama, K. Kishio, H. 
Ikeda, R. Yoshizaki, Z. Hussain, and Z.-X. Shen, Phys. 
Rev. Lett. 85, 2581 (2000). 

We want to emphasize that the dispersion is indeed badly 
known because of the low resolution of ARPES measure- 
ments. Even in the long studied compound BSCCO, it 
was until recently that some kink in the dispersion is re- 
alised. The kink results in a reduction of velocity by a 
factor of 2 or more. It is interesting to note that BSCCO 
is a compound structurally similar to YBCO, and we ac- 
tually need a vf of YBCO several times smaller than that 
of LSCO. 

J.P. Lu, Phys. Rev. Lett. 68, 125 (1992). 

A. Abrikosov, Phys. Rev. B 62, 15156 (2000). 

T.E. Mason, G. Aeppli, and H.A. Mook, Phys. Rev. Lett. 

68, 1414 (1992). 

A. A. Abrikosov, Phys. Rev. B 57, 8656 (1998). 

V.M. Krasnov, A. Yurgens, D. Winkler, P. Delsing, and 

T. Claeson, Phys. Rev. Lett. 84, 5860 (2000). 

T. Sato, T. Yokoya, Y. Naitoh, T. Takahashi, K. Yamada, 

and Y. Endoh, Phys. Rev. Lett. 83, 2254 (1999). 

K. Gorny, O.M. Vyaselev, J.A. Martindale, V.A. Nandor, 

C.H. Pennington, P.C. Hammel, W.L. Hults, J.L. Smith, 

P.L. Kuhns, A.P. Reyes, and W.G. Moulton, Phys. Rev. 

Lett. 82, 177 (1999). 

K.-K. Voo, H.-Y. Chen, and W. C. Wu (unpublished). 
J. Rendell and J. Carbotte, Phys. Rev. B 53, 5889 (1996). 
K.A. Kouznetsov, A.G. Sun, B. Chen, A.S. Katz, S.R. 
Bahcall, J. Clarke, R.C. Dynes, D.A. Gajewski, S.H. Han, 
M.B. Maple, J. Giapintzukis, J.-T. Kim, and D.M. Gins- 
berg, Phys. Rev. Lett. 79, 3050 (1997). 
R. Kleiner, A.S. Katz, A.G. Sun, R. Summer, D.A. 
Gajewski, S.H. Han, S.I. Woods, E. Dantsker, B. Chen, 
K. Char, M.B. Maple, R.C. Dynes, and J. Clarke, Phys. 
Rev. Lett. 76, 2161 (1996). 

A.G. Sun, D.A. Gajewski, M.B. Maple, and R.C. Dynes, 
Phys. Rev. Lett. 72, 2267 (1994). 

S. Quinlan and D. Scalapino, Phys. Rev. B 51, 497 
(1995). 

L. Yin, S. Chakravarty, and P. W. Anderson, Phys. Rev. 
Lett. 78, 3559 (1997). 

Very recently there is a report of a compreh ensive 



ARPES meas urement on LSCO by A. Ino et al. in cond 



mat/0005370 



As debates on ARPES measurements fre- 
quently occur, such as the FS is hole- or electron-like in 
BSCCO [see P.V. Bogdanov et al 



:ond-mat/0005394 



and references therein], and the existence or not of an 
extended-vHS in Sr2Ru04 [see A.V. Puchkov et al., Phys. 
Rev. B 58, R13322 (1998), and references therein], we 
think that it is still early to discuss the effect of the re- 
sult by A. Ino et al. on our conclusion. 

A.V. Chubukov, B. Janko, and O. Tchernyshyov, cond 



mat/0012065 



D.S. Marshall, D.S. Dessau, A.G. Loeser, C-H. Park, 
AY. Matsuura, J.N. Eckstein, I. Bozovic, P. Fournier, 
A. Kapitulnik, W.E. Spicer, and Z.-X. Shen, Phys. Rev. 
Lett. 76, 4841 (1996). 



13 



